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ABSTRACT 

The problem of a non-steady planar radiation mediated shock (RMS) breaking out from a surface with a 
power-law density profile, p cx x", is numerically solved in the approximation of diffusion with constant opacity. 
For an appropriate choice of time, length and energy scales, determined by the breakout opacity, velocity and 
density, the solution is universal, i.e. depends only on the density power law index n. The resulting luminosity 
depends weakly on the value of n. An approximate analytic solution, based on the self-similar hydrodynamic 
solutions (jSakurai 1960) and on the steady RMS solutions (e.g. Weaver 1976), is constructed and shown to 
agree with the numerical solutions as long as the shock is far from the surface, t ^ c/\s,-n. Approximate analytic 
expressions, calibrated based on the exact solutions, are provided, that describe the escaping luminosity as a 
function of time. These results can be used to calculate the bolometric properties of the bursts of radiation 
produced during supernova (SN) shock breakouts. For completeness, we also use the exact breakout solutions 
to provide an analytic approximation for the maximum surface temperature for fast (Vsh ^ 0.1) non-thermal 
breakouts, and show that it is few times smaller than inferred based on steady-state RMS solutions. 

Subject headings: radiation hydrodynamics — shock waves — supemovae: general 



1. INTRODUCTION 

During a SN explosion, a strong shock wave propagates 
through and ejects the stellar envelope. This shock moves 
with high velocities, /?sh = Vsh/c > 0.01, through plasma with 
densities p ^ 10~'° - 10~^g cm""'. Under such conditions the 
energy density is dominated by radiation, and the opacity is 
dominated by electron scattering. The kinetic energy is con- 
verted to thermal energy through Compton scattering of free 
electrons. The transition layer of such a Radiation Medi- 
ated Shock (RMS) has a finite width, with optical depth of 
5t ^ c/vsh (Weaver 1976). 

As long as the shock is deep within the envelope, where 
the density does not vary significantly over the shock width, 
its propagation is well described by the (ideal fluid) hydro- 
dynamics equations, considering the shock as a discontinuity. 
For envelope density that varies as a power-law of the dis- 
tance from the stellar edge, as is the case for a wide range of 
progenitors, the hy drodynamic pro files are given by the self 
similar solutions of ISakurail (I1960I) . Once the shock reaches 
a distance from the surface which is comparable to the width 
of the front, i.e. an optical depth t from the surface for which 
T ^ c/vsh, the hydrodynamic solutions are no longer valid. At 
this point, the radiation in the shock transition region diffuses 
out a n d a burst of radiation is e r nitted (Colgate 1974; Falk 
19781; iKieiiT&Chevaiie? 119781; 'Ensman & Burr owa ■1992i; 
Matzner & McKee 1999; Blinnikov et al. 2000; iKatz et alj 
2009c iPiro et al.ii2010c iNakar & Sariii2010i) . The light curve 
and spectrum of the escaping radiation depends on the struc- 
ture of the shock transition. 

The structure of steady state RMSs propagating through 
a homogeneous medium w ith constant ve locity was solved 
both in the non-relativistic (IWeaveiill976l) and in the highly 
relativistic dBudnik et al.l 120101) regimes. An approximate 
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analytic descript ion of the shock structure was derived 
jKatz et alj|2009t) for fast shocks, /3sh > 0.1, where the radi- 
ation is in C ompton equilibrium and far from thermal equi- 
librium (Weaver 1976). These solutions do not provide, how- 
ever, an accurate description of the shock near breakout, since 
at breakout the density varies on a scale comparable to the 
shock width and the shock begins to lose energy to the out- 
flow of radiation. 

The emitted burst was previously studied either by nu- 
merically calculating the early emission from specifi c pro- 
genitors (includ ing Ensman & Burrows 1992; Blinniko v et al.l 
l200^, 'Utrobin' '" 2007i; ,Tominaga et al.. ,2009c ,Tolstovl l2010r 
Dessart & Hillier n2'010l) . or by providing analytic order of 
magnitude estimates for the breakout properties based on 
the ste ady state solutions (e.g. Matzner & McKee 199%, 
IKatz et al. 2009; Nakar & Sari 2010). The bolometric proper- 
ties of a RMS breaking out from a n envelope with an expo nen- 
tial density profile was solved bv lLasher & Chaiil ([T97l (We 
note that the reported results were restricted to e arly times, 
in which only 0.5 of the energy is emitted, see § 14 3b . The 
accurate time dependent spectrum of shock breakout from a 
general progenitor has not yet been calculated. 

For non-relativistic breakouts, the transport of radiation is 
well described by diffusion. As argued below, in § |4] the dif- 
fusion approximation holds at least for Vsh/c < 0.3. More- 
over, for such velocities the maximum electron temperature 
reached is ^ nieC^ (see fig. |9]l, which implies that the scat- 
tering cross section is frequency independent. Thus, the hy- 
drodynamic and radiation energy density profiles may be ob- 
tained by solving the hydrodynamics equations coupled to a 
constant opacity diffusion equation for the radiation energy 
density (iLasher & Chan 1979). In this paper, the time depen- 
dent hydrodynamic and radiation flux profiles before, during 
and following breakout are derived by solving these equations 
(neglecting plasma pressure). 

The exact formulation of the problem is given in § |2] the 
numerical solutions are described in § |3] and the results for 
the hydrodynamic profiles and emitted flux and energy (per 
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unit area) are presented in § B] These results are used in the 
second paper of this series (iKatz et al.ll201 Ih to find the hy- 
drodynamic and bolometric properties of SN shock breakouts 
taking into account the finite radius of the progenitor and the 
lowest order relativistic corrections. In § 2] we explain how 
the exact solutions obtained here may be used to derive the 
spectral properties of the breakout. The temperature profiles 
and the spectral properti es of the br eakout are described in a 
third paper in this series (ISapir et al. 2011). For completeness, 
we give in §|5]the maximum surface temperature obtained for 
various breakout velocities and densities. Our main conclu- 
sions are summarized in § |6l 

2. FORMULATION OF THE PROBLEM 

Consider the problem of a planar RMS propagating in the 
positive X direction through a plasma occupying the a; < re- 
gion. The initial, t = -oo, density of the plasma decreases 
towards the initial position of the surface, x = 0, as 



p|/=-oo OC 



\x\" x<0 
x>0 



«/(«+!) 



(1) 



T(x,t) is the optical depth of the matter lying ahead of x, 

/>00 /"OO 

T = npdx = K pdx, (2) 



where k is the constant opacity. 

The analysis is preformed under the following assump- 
tions/approximations : 

• The velocities are non relativistic, /3 = v/c <C 1; 

• The internal energy (and pressure) of the matter are ne- 
glected; 

• Photon transport is described by diffusion with a con- 
stant opacity k. 

Asymptotically, at large optical depth from the surface, the 
hydrod ynamic profiles are described by the self-similar solu- 
tions of lSakurail (I1960I) . The shock velocity follows a power- 
law dependence on the optical depth Vsh oc r"'^^" '''^/*"^'^ with 
A depending on the density profile index n and on the value 
of the adiabatic index 7. In the case of radiation domination, 
7 = 4/3 and A(n) « 0. 19n (Grasberg 198 1). In the limit n^oo 
the problem co rresponds to a strong shock in an exponential 
density profile (: Haveslll968l) . 

The shock has a finite width St ^ /3~^l (e.g. IWeaveijri976l) 
where /3sh = Vsh/c. At an optical depth of r ^ f3~^, energy 
diffuses through the surface and the shock dissolves. Note, 
that at any finite optical depth, the hydrodynamic profile is 
"smoothed", i.e. the shock transition layer is not a discontinu- 
ity, due to the diffusion of energy, and the "shock velocity" is 
not well defined. 

2.1. Hydrodynamic and energy diffusion equations 

The following Lagrangian coordinates are used to describe 
the flow. For each mass element, m = -t/k is minus the mass 
per unit area enclosed between the element and the surface 
while x,\ and e are the position, velocity and energy density 
of the element respectively. In this problem, all mass ele- 
ments start with negative positions and throughout the evolu- 
tion move outward with positive velocities. 



Under the above assumptions, the equations for x, v and e 
as a function of m and the time t are 



d,x = V, 

d,W = -dmp, 

d,(e/p) = -d,„j-pd,„v, 



J 



3k 
e = 3p 



(3) 
(4) 
(5) 

(6) 

(7) 



where p = d^m is the density, j is the energy flux (energy cur- 
rent density), p is the radiation pressure and k is the constant 
opacity. 

The initial density profile and the asymptotic shock velocity 
at large optical depth can be parameterized by 



and 



Pm = P0(/30T)'"'("^') 



Ah ^ A)(A)T)-^/("^'^ 



(8) 



(9) 



respectively. Note that the equations are normalized so that 
the density and shock velocity at t = /?g ' woul d equal pn and 
A) respectively if the hydrodynamic solution (ISakurailll960h 
was applicable to this point. That is, breakout occurs at t ~ 

At' with Ah ^ A- 

Radiation is assumed to freely escape from the surface. 
Given that the smallest optical depth over which changes can 
occur is much larger than unity, r ^ /?g ' ^ 1, the energy den- 
sity e is set to zero at the surface, 

e(A)T = 0) = 0. (10) 

The velocity at the surface is determined by equations (IS)-© 
as follows. Combining equations Q,® and (U the following 
general relation is obtained (ILasher & Chanlll979lF 1 



d,v= 

c 



(11) 



Since the energy flux is continuous across the surface, the ve- 
locity of the surface is related to the flux by (ILasher & ChanI 
Il979h 



5,v(At = 0)=-7(At = 0). 



(12) 



The dimensionful parameters k, po and vo = A)C along with 
the dimensionless parameter n cornpletely define the problem 
through the initial density profile (|8]l, the asymptotic velocity 
at early times (|9|l and the boundary conditions (fTOt and (fT2l i. 
along with the equations of motion ©-(iTli. 

2.2. Dimensionless equations 
By introducing dimensionless variables, A=A /Aq with 
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1 



Xq = 



'"0 = — 5-' 



to=—, po = PoVo, (13) 

Vo 



Note that Eq. {TT) is a direct consequence of tlie assumption that the 
plasma's pressure is negligible. This equation is equivalent to stating that in 
scattering events of forward/backward symmetiy, the momentum of the radia- 
tion is completely transfen'ed to the mediu m. This r elation holds independent 
of the angular distribution of the radiation dLasher & Chan 1979). 



the following dimensionless equations are obtained: 

V =djx, 
dfw=-dmp, 

with an initial density profile 

An = l-r/^"'^ 

asymptotic shock velocity 

I - |-A/(n+l) 
Vsh >\m\ " 

m— >-oo 





^^peak/^peak^ ^^peak/^peak^ 




box size /l" 
box size X ^ 


1.1x10"-^ 1.2x10"-^ 
1 x 10-3 5 X 10-3 


4.8 X 10-3 2.3 X 10-3 
4x10-* 1.2x10-2 


resolution /tF 
resolution 


1.4x10-^ 2.8x10-^ 
3 x 10-3 1.3 X 10-2 


2.3 X 10^ 8.6 X 10"^ 
7 X 10-5 5 5 y, iQ-A 



(14) 
(15) 
(16) 



and surface boundary conditions, 

p(0) = 0, dTm = -d^p\M=o- (17) 

The system of equations (fl4]l-(fT7]i has one free parameter, 
namely the value of the initial density power-law index n. 
Therefore, for a given value of n, the time dependent profiles 
of the dimensionless quantities \,p,x,p are universal (inde- 
pendent of Vq, pq and k). This universality is not surprising. 
For a given value of n, the (dimensionful) problem is defined 
by the values of three dimensionful parameters, namely vo, po 
and K as well as the speed of light c. The only dimensionless 
parameter (aside from n) that can be obtained is v/c, which is 
taken to zero in the non relativistic limit. 

3. NUMERICAL INTEGRATION 

Equations ( fT4l i were numerically solved using an implicit 
scheme in time, a grid with uniform spacing in mass, and 
using constant time steps. Any boundary condition on the 
inner boundary, combined with an initial power law density 
profile would result in an accelerating shock which would ap- 
proach the Sakurai solution, (fTST l. far from the boundary, as 
long as the box is large enough. In order to achieve fast con- 
vergence with box size, the initial distributions of positions, 
velocities and pressure were chosen to be an approximation 
of the solution at an early stage when the shock was close to 
the inner boundary. The precise profiles were described by 
an ansatz combining the steady state structure of a radiation 
mediated shock with the time dependent Sakurai self-similar 
solution (see § |A]i. The velocity and energy flux of the inner 
boundary (as a function of time) were set to be equal to the 
values given by the Sakurai solution. The resulting profiles 
converged quickly with box size as shown below. 

As a consistency check, a set of different boundary condi- 
tions was used in an additional calculation. In this case, the 
inner boundary moves with a constant velocity and energy is 
deposited at the innermost cell at the initial time. The result- 
ing profiles were found to converge to the same profiles that 
were obtained with the Sakurai-Weaver ansatz, albeit much 
slower 

The hydrodynamic profiles shown below were obtained 
with a 6000 cells box of width dm = 0.08 and time steps of 
c/f w 5 X 10"^. Initially, the shock was positioned at a (mass) 
distance of (5m w 10 from the inner boundary. The emitted en- 
ergy flux and fluence results were obtained with a 1200 cells 
box of width dm = 0.025 and time steps of t/f ss 1 x lO"'*. Ta- 
ble [Upresents the convergence with increasing resolution and 
box size of paramet ers d escribing the characteristics of the 
flux and fluence (see l4.3l for definitions). As indicated by the 
table, numerical results presented for the values of these pa- 
rameters are converged to an accuracy of better then 1%. The 



"See table[2]for definition and resulting value 

''With respect to the nominal calculation described in the text 



TABLE 1 

Convergence of the numerical solution 



table presents convergence with respect to box size and reso- 
lution. The convergence with respect to changes in time step 
size is better than 10"^ and not presented. 

4. RESULTS 

The resulting profiles of the density, velocity and pressure 
as a function of normalized mass at different times are given 
in § 14.11 for n = 3 . The shock velocity at different times is 
calculated using different prescriptions, a nd c ompared to the 
hydrodynamic self-similar velocity, in § 14.21 The resulting 
emission of radiation from the surface and the asymptotic dis- 
tribution of velocities is presented in § 14.31 t /to = is chosen 
as the time when the Sakurai shock reaches the surface, ac- 
cording to the integration of (|9}. 

4. 1 . Density, velocity and pressure profiles 

The hydrodynamic profiles of the density, pressure and ve- 
locity of the gas as the shock progresses to the surface are 
shown in figures [T] |2] and [3] respectively. The results of 
the simulation are compared to a Sakurai-Weaver ansatz that 
combines the Sakurai accelerating shock with the Weaver 
steady state RMS, without any fitted parameters (see Ap- 
pendix). As can be seen, the ansatz agrees with the results 
to a high accuracy up to the time when the shock reaches an 
optical depth of t ^ /?g ' . 



10' 



10' 



10" 




Fig. 1. — The noiTnalized density profiles as a function of the normalized 
mass (optical depth) at different times prior to breakout (blue solid lines). 
Also plotted are the density profiles of the SW ansatz at the appropriate times 
(red dashed lines, see § |A). The arrow indicated the shock propagation di- 
rection. I3qc is the shock velocity "at breakout", i.e. t he v elocity for which 
r = 1//3 in the Sakurai solution, and to is defined in eg. 1131 The peak emitted 
flux is obtained at ?pcak/'o = —1-25. 

The profiles at different times after breakout are plotted in 
figures |4]|6] The results of the simulation are compared to the 
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Fig. 2. — The normalized pressure profiles as a function of the normalized 
mass at different times prior to breakout. Line styles/colors are as in figure[T] 




10° 10' 10' 

Fig. 3. — The normalized velocity profiles as a function of the normalized 
mass at different times prior to breakout. Line styles/colors are as in figure[T] 

Sakurai solution of the pure hydrodynamic equations. Note 
that the SW ansatz does not apply at these times, which are 
greater than the breakout time, since a shock no longer exists. 
As can be seen in the figures, diffusion affects the profiles only 
close to the boundary r < /3q'. This is expected since the dif- 
fusion depth has a const ant optical depth Tdiff ~ through- 
out the expansion (e.g. [ Matzner & McK ee 19991; iPiro et all 
l2010l:lNakar & Sarill2010l) . The difference between the Saku- 
rai solution and the numerical solution at the point t = /3q' is 
^ 10%. Note that in the case « = 0, the optical depth of the 
diffusion depth grows with time and the deviations from the 
Sakurai solution are observed on this growing scale. 

The results of the velocity profiles also show for which ve- 
locities the diffusion approximation is appropriate in describ- 
ing the radiation transport. The diffusion approximation is 
valid as long as the velocity does not change significantly 
over a photon mean free path, and this requirement can be ex- 
pressed as drP <C 1. During propagation, the velocity deriva- 
tive increases as the shock approaches the surface and the 
maximal value is dr(3 ^ 1.0/3(). Therefore, we conclude that 
the diffusion approximation is likely valid for /3o ^0.3. 




Fig. 4. — The normalized density profiles as a function of the normalized 
mass at different times following breakout (blue solid lines). Also plotted are 
the density profiles of the self-similar hydrodynamic Sakurai solution at the 
appropriate times (black dashed lines). 
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Fig. 5. — Normalized pressure as a function of the normalized mass at 
different times following breakout. Line styles/colors are as in figure|4] 




Fig. 6. — Normalized velocity as a function of the normalized mass at dif- 
ferent times following breakout. Line styles/colors are as in figure|4] 
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4.2. Shock Velocity 

Strictly speaking, the position and the velocity of the shock 
are not well defined, since the diffusion of radiation smooths 
the hydrodynamic shock discontinuity. Here, different pre- 
scriptions for defining the shock velocity given the hydro- 
dynamic profiles are used for comparing the results to the 
pure hydrodynamic Sakurai solution (which has a well de- 
fined shock position and velocity). The effective shock po- 
sition msh,eff(f) is defined as the point where the density in- 
creased by some chosen factor 1 < / < 7 compared to the 
original density, 

/0(msh.eff , f ) = fPin (wJsh.eff) ■ (18) 

In this section, / is chosen to be / = \/7. The effective shock 
velocity is obtained by the following two prescriptions which 
reduce to the correct shock velocity in the steady state case. 

1 . The effective shock velocity is set to be the time deriva- 
tive of the effective shock position, 

vl,ff(f)= -^x(m,h,eff,0- (19) 




Fig. 7. — The normalized effective shock velocities as a function of the 
effective position of the shock (in terms of nomialized mass). The blue and 
red solid Unes are the effective shock velocities estimated from equations | |19I 
and 120 1 respectively and the black dashed line is the Sakurai solution. The 
separation of the red and blue curves is an indication of the deviation from 
the (ideal fluid) hydrodynamic evolution. 



2. Based on the steady state shock relation in the up- 
stream frame, p = PupVVsh (which holds for a steady 
state shock moving with velocity Vsh through a homo- 
geneous medium with pup), the effective shock velocity 
is set to be 



Vsh.eff(0 = 



Pin('«sh,eff)v(OTsh,eff,0' 



(20) 



The numerical shock velocity is compared in figure|7]to that 
obtained in the self-similar (ideal fluid) hydrodynamic solu- 
tion. As can be seen in the figure, when the shock is far from 
the surface, the two definitions of the effective shock velocity 
agree with each other and with the Sakurai solution. The two 
definitions lead to different estimates of shock velocity close 
to the surface. This is an indication of the deviation from the 
(ideal fluid) hydrodynamic evolution, which is expected to oc- 
cur when the shock width is comparable to the distance from 
the surface. 

4.3. Emission of radiation from the surface 

We next consider the energy flux (i.e. luminosity per unit 
area) emitted from the surface. Lit) = jifioT = 0,f), and the 
amount of energy per unit area emitted up to a given time f 



8{t)= / C{t')dt\ 



(21) 



during and following the shock breakout. The emitted flux 
and energy, normalized to 



(22) 



are shown in figure |8] As can be seen in the figure, the emit- 
ted flux rises sharply and later decreases, asymptotically ap- 
proaching a power law C{t) oc f"^/^. The amplitude and posi- 
tion of the peak and the full width at half maximum are given 
in table 12] 

The late time power law decline, £(f) oc t~^l^, is expected 
for any density profile index « > dMatzner & McKeelll999t : 
Firo et al...2010;.Nakar & Sari2010r) . To see fliis note fliat the 



diffusion optical depth is constant in time and the luminos- 
ity is proportional to the pressure at r '-^ /3"^ which roughly 
decreases adiabatically p cx p^l^ with the density decreasing 
linearly with time. Note that in the case « = the diffusion 
optical depth grows with time. Taking this into account, the 
asymptotic emitted flux follows C{t) cx f"'^-"^ for this case. The 
dimensionless emitted energy and energy flux can be approx- 
imated at late time (f > fpeak + flf^o) by 



t ^peak 
a,t() 



-1/3- 



• (f > fpeak + fl/fo), 



^oo f ^ ^peak 



3aft{) \ Gito 



-4/3 



• (t > tpeak + a,to), 



(23) 



where fpeak is the time at which C(t) reaches its maximum, a, 
is a free parameter, and (x > y) is a function that is unity for 
X > y and zero otherwise. The fitted values of a, and £00 for 
n = 3/2, 3 are given in table |2] The corresponding power-law 
fit for the case « = 3 is shown in figure [8] 

The early stages of the rise in emitted energy flux are due to 
photons diffusing ahead of the shock while it is still far from 
the surface. Since the distribution of these photons drops like 

cx e"^ away from the shock, the rise in the emitted energy flux 
can be approximately described by 



£(0 = 



(24) 



The fitted values of £i,fli and b; are given in table |2] and the 
corresponding fit for the case n = 3 is shown in figure [8] 

Using equation ( fT2b . the acceleration and velocity of the 
surface (outermost mass element) are related to the emitted 
energy flux and fluence by 



and 



dM/3oT = Q,t)=-m 
c 



wWoT = 0,t) = -£(t) 
c 



(25) 
(26) 
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Peak characteristics 


Eq. (23) Fit 


Eq. (24) Fit 


n" 




£oQ^£o = Vcx) /vo 0/^ 




3 

3/2 


0.72 -1.25 1.61 
0.77 -0.78 1.49 


2.03 0.1 
2.13 0.129 


1.25 1.11 1.57 
0.45 1.08 1.14 



"Density power law index, Eq. {T) 
''Peak Luminosity 

'^Time of peak luminosity, measured from time of expected breakout 
Luminos ity full width half maximum 
■^See Eq. (23j 
'See Eq. (24) 

TABLE 2 

Light curve characteristic values 



respectively. Thus the plotted normalized emitted flux and 
fluence are, respectively, also the normalized acceleration and 
velocity of the surface 

9,v(/3oT = 0,f)ro/vo = /:(f)//:o, 
v(A)r = 0,f)/vo = £(f)/fo, (27) 
and the surface velocity can be approximated by 

-1/3-1 



v(/3oT = 0,f): 



1- 



t - fref 
fl,f() 



•(f>W) (28) 



where Voo/vo = Soo/^o is given in table |2] 

For comparison, the velocity of the mass element at r = /3q ' 
in the pure hydrodynamic Sakurai solution is also plotted in 
figure[8](dashed dotted lines. Note that the velocity of the sur- 
face in the Sakurai solution is for f < and oo at f > 0). 
As can be seen, the pure hydro solution can be used to ap- 
proximate the emitted radiation and the surface motion at late 
times. In particular, the asymptotic velocity achieved by the 
surface is equal to the maximal velocity achieved in the Saku- 
rai solution (for the mass element at r = /3q') to an accuracy of 
^ 1%. This implie s that previ ous estimates of the velocity of 
the fastest element (iMatzner & McKee 1999) are accurate. As 
far as we can tell, this high level of accuracy is coincidental. 




Fig. 8. — Nomialized (see eg. 1221 emitted energy flux (solid blue line) and 
energy per unit area (solid red line) as a function of time since (expected) 
breakout. These are equal to the noiTnalized acce leration and velocity of the 
outermost mass element (see discu ssion in § |4.3| . Also plotted are the fits to 
the emitted flux given by Eq. (23) for t/to> and by Eq. (24) fort /to < -1 
(magenta dashed-dotted lines) and the velocity of the r/3o = 1 mass element in 
the Sakurai solution (dashed black line). The inset shows the same quantities 
on a larger time scale. 



5. TEMPERATURE 

In a following paper (I Sapir et aDl201 Ih . the spectral com- 
ponents of the breakout burst are calculated by extending the 
planar calculation to include a determination of the temper- 
ature. This is done by solving a second diffusion equation 
for the photon number density. Assuming Compton equilib- 
rium, the temperature can be obtained from the pressure and 
photon number density. In particular, this method allows the 
first accurate calculation of the spectra for fast shocks [3 > 0.1 
where the radia t ion is far from equ ilibrium (IWeaverlll976t 
iKatz et alJl2009l:lNakar & Sarill2010l) . This can be compared 
with the analytical , order of magnitude estimates of peak tem- 
perature (iKatz et al. 2009; Nakar & Sari 2010), and expected 
spectral features (.Nakar & Sari,,2010,) . For illustration, the re- 
lation between peak temperature T^ax, Po and Vo is shown in 
figure|9]for breakout in an envelope consisting of pure hydro- 
gen or pure helium (the results depend on the atomic number 
and charge through the combin ation 7? jA) . Also plotted is 
the prediction of equation (18) in lKatz et alj ( 120091) for steady 
state RMS. The resulting temperature deviates from the value 
achieved by a steady state shock propagating in a density po 
with velocit y vp, used as an approximation for the break out 
temperature (IKatz et alJl2009l;]Nakar & Sarill20Tol) by factors 
of a few. 

The relation between /3o, po and Tysko. can be fitted by, 

/3i = alogio(rkev) + ^logio(7kev) + clogio(p_9) + J (29) 

where T = Tkev keV and p = 10~^p-9 g cm"^. For « = 3, rele- 
vant to BSGs and WRs where the velocities can reach > 0. Ic, 
the best fit parameters are a = 0.03, /? = 0.1 33 , c = -0.0267 and 
d = 0.153. The photon energies at the spectral peak are ap- 
proximately hiy « 3rtnax- The detailed pro perties of the spec - 
trum of the emitted burst will be studied in (ISapir et al.ll201 Ih . 

6. DISCUSSION 

The problem of a planar RMS breaking out of a gas with 
a power-law density profile and with constant opacity was 
solved in the diffusion approximation. The main results are 
summarized below. 

• For a given density power law index n the emitted en- 
ergy flux and the hydrodynamic behavior have univer- 
sal profiles in space and time, up to scaling of density, 
velocity and opacity. The scaling of time, distance, 
mass and energy density are given in Eq. ([13) (in this 
equation, fioc is the shock velocity "at breakout", i.e. 
the velocity for which r = 1 / /3sh in the Sakurai solution, 
and Po is the initial density at r = 1 //3o). 

• The emitted energy flux can be fitted by a super expo- 
nential rise, jO(xe~' , at early times and by a power law 
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Fig. 9. — The shock velocity i3q as a function of maximum surface tem- 
perature for different values of breakout density po. Crosses represent the 
calculated value, solid cur ves show the curv ed fit and dashed curves show the 
results of equation (18) of lKatz et aH 120091) for steady-state RMS. 

decay, C cx t~^^^, at late times, see Eqs. ( l24b and ( l23T l 
(for n = the late time luminosity follows C oc f"''/**). 
The values of the parameters that appear in these eqs., 
along with properties of the emitted flux peak, are pro- 
vided in table |2] (the normalization is given in eq. l22T l. 
The emitted energy flux as a function of time for power 
law indexes in the range n = 1 - 10 is fully provided in 
tableland shown in figure [TOl 

• The luminosity depends weakly on the decreasing den- 
sity structure (e.g. on the value of «). For power law 
indexes in the range n = 1-10 the luminosity changes 
by less than 25% (see table[3]and figure [TO]i. 

• The acceleration (velocity) of the different mass ele- 
ments is linearly related to the energy flux (fluence) 
passing throu gh the mass element b y d,v = k/cJ (see 
Eq. (in), and lLasher & Chanlll979l) . In particular the 



velocity of the surface, which is also the fastest mov- 
ing element, is given by the fluence emitted from the 
surface. 

• At late times, t > 2to, the emitted energy flux and flu- 
ence are approximately equal to the acceleration and 
velocity at r = c/vo in the Sakurai solution (see figure 

B- 

• The hydrodynamic profiles before breakout t < -3fo 
(see figures [T] |2] and O can be accurately described 
by an ansatz (see § |A]i combining the hydrodynamic 
power-law solutions of Sakurai and RMS steady state 
solutions of Weaver. 

• There is no density discontinuity at late stages, suggest- 
ing that a coUisionless shock does not form in a planar 
breakout. The formation of the dense shell observed in 
(lEnsman & Burrowslll992l) . if real, must be related to 
other physical processes (e.g. spherical expansion or 
transport of radiation at r < 1). 

• We explained in §|5]how the exact solutions obtained 
here may be used to derive the spectral properties of 
the breakout. The temperature profiles and the spec- 
tral properties of t he breakout are described in detail in 
(ISapir et al.ll201 ll) . For completeness, we gave in § |5] 
the maximum surface temperature obtained for various 
breakout velocities and densities, see fig.|9]and eq. ( |29] l. 
For fast (Vsh 0. 1) non-thermal breakouts the tempera- 
ture is few times smaller than inferred based on steady- 
state RMS solutions. 
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APPENDIX 

PROFILE ANSATZ COMBINING A TIME DEPENDENT HYDRODYNAMIC SOLUTION AND THE STEADY 

STATE RMS SOLUTION 

In this section we describe an ansatz that approximates the profile of an RMS propagating through a cold non homogenous 
medium with initial density profile /3in(m). We found this ansatz to accurately describe the exact profiles for the case of an initial 
power law density distribution (see figures [T]|3). 

It is assumed that the solution of the pure hydrodynamic problem of an infinitely thin, strong RMS propagating through the 
medium is known. Velocities are measured in the laboratory frame, where the unperturbed medium is at rest. The hydrodynamic 
solution is described by the shock position and velocity as a function of time, msh(f) and Vsh respectively, and by the values of 
the position, pressure and velocity of every mass element m, XH(ni,t),pH{m,t),VH(m,t) respectively. The density is given by 
PH(ni,t) = dx^m. Denote the unperturbed density immediately upstream of the shock by Psh(^) (Psh = Pin('Msh(0)- By assumption, 
for m > nish we have pnim^t) = VH{in,t) = 0. The ansatz is given by 

JCanz('«, t) = XH + (K/9sh/3sh)~'^('«), 



,(m,f) = < 6 -X 

\ 7VshV(m) m > m^h 

PH- v(m) m < nish 



where 



m= K/3sh(m-msh) 



(Al) 



(A2) 
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and i, V are given by 

- In m ■ (m > 0), 

7 V 2 y 7 ' 

1 

(A3) 

Eq. ( lAU is motivated by the fact that the spatial width of the shock is Sx ^ (Kpsh/3sh)~' and the mass within the shock region 
is Sm ^ (K/3sh)~'- The expressions for x and v where uniquely determined by the following requirements: 1. for a steady state 
shock propagating through a homogenous medium, this ansatz should exactly reduce to the known analytical solution, 2. at the 
shock and far from the shock the spatial deformation goes to zero 

X > 0. (A4) 

(5m— ^itoo.O 

A comparison of this ansatz to the exact numerical results for the case of an initial density distribution p cx jmp^^ (« = 3) is 
presented in figures [T]|3] 

NUMERICAL VALUES OF £(r). 

The numerical results for the instantaneous luminosity per unit area, C{t), resulting from the numerical planar solution, are 
given in tableland are plotted in figure [TO]. Time is measured with respect to the time of peak luminosity. The differences in 
luminosity between the « = 1-10 cases and the « = 3 case are less than about 25% at the time interval -1 < (f-fpeak)/fo < 20, 
where about 80% of the energy is emitted. 



(;-;peak)A(il 


C/C(^ (n = 3t 


£(n= l)/£(n = 3) 


/;(« = 3/2)/£(n = 3) 


C{n = 10)/Z:(n = 3) 


-2.00 


0.00 ife 


1.45 


1.26 


1.26 


-1.75 


0.00712 


1.14 


1.08 


1.34 


-1.50 


0.0236 


0.979 


0.984 


1.34 


-1.25 


0.0665 


0.908 


0.941 


1.28 


-1.00 


0.157 


0.903 


0.939 


1.17 


-0.75 


0.31 


0.949 


0.968 


1.06 


-0.50 


0.5 


1.02 


1.01 


0.962 


-0.25 


0.66 


1.1 


1.05 


0.902 


0.00 


0.719 


1.13 


1.07 


0.882 


0.25 


0.669 


1.1 


1.06 


0.899 


0.50 


0.557 


1.04 


1.02 


0.939 


0.75 


0.435 


0.99 


0.993 


0.989 


1.00 


0.333 


0.957 


0.974 


1.03 


1.25 


0.257 


0.946 


0.967 


1.06 


1.50 


0.201 


0.948 


0.969 


1.07 


1.75 


0.162 


0.958 


0.977 


1.07 


2.00 


0.133 


0.97 


0.986 


1.06 


2.25 


0.112 


0.983 


0.997 


1.05 


2.50 


0.0963 


0.996 


1.01 


1.04 


2.75 


0.0839 


1.01 


1.02 


1.02 


3.00 


0.074 


1.02 


1.02 


1.01 


3.25 


0.0661 


1.03 


1.03 


1.0 


3.50 


0.0595 


1.04 


1.04 


0.99 


3.75 


0.0541 


1.04 


1.04 


0.981 


4.00 


0.0494 


1.05 


1.05 


0.973 


4.25 


0.0455 


1.06 


1.05 


0.966 


4.50 


0.042 


1.06 


1.06 


0.96 


4.75 


0.039 


1.07 


1.06 


0.954 


5.00 


0.0364 


1.08 


1.06 


0.95 


6.00 


0.0285 


1.1 


1.08 


0.934 


8.00 


0.0194 


1.14 


1.09 


0.913 


10.00 


0.0144 


1.17 


1.11 


0.9 


15.00 


0.00845 


1.22 


1.13 


0.879 


20.00 


0.00578 


1.26 


1.14 


0.867 


30.00 


0.00338 


1.3 


1.16 


0.851 


50.00 


0.00172 


1.33 


1.18 


0.836 


75.00 


0.001 


1.36 


1.2 


0.825 


100.00 


0.000687 


1.38 


1.21 


0.817 


125.00 


0.000513 


1.41 


1.22 


0.811 



^Time relative to t^^-^i^, the time when C(t) peaks, nomialized to tg = c/(kpqVq). 

''Energy flux emitted from the surface, normalized to Co = PoVq (For a surface of 4itR^ this is equivalent to the luminosity normalized to 4itR^ po\^). 
'^Density is assumed to increase away from the surface as a power law p cxx". 
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Fig. 10. — Normalized emitted energy flux as a function of the normalized time relative to time of peak emitted energy flux, for different values of the power-law 
index n. Also plotted is the fitting function C[t) = 0.33£o [(f~'pcak) Z'o]^^^ (dashed black line) corresponding to l |23t . 
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